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Non-relativistic bound states at finite temperature (II): the muonic hydrogen 
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We illustrate how to apply modern effective field theory techniques and dimensional regularization 
to factorise the various scales which appear in QED bound states at finite temperature. We focus 
here on the muonic hydrogen atom. Vacuum polarization effects make the physics of this atom at 
finite temperature very close to that of heavy quarkonium states. We comment on the implications 
of our results for these states in the quark gluon plasma. In particular, we estimate the effects of a 
finite charm quark mass in the dissociation temperature of bottomonium. 
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(SI . I- INTRODUCTION 

In a previous paper [l| we showed how to apply modern effective field theory (EFT) techniques to the hydrogen 
atom at finite temperature. They provide a systematic way to separate the physics occurring at the various dynamical 
scales involved in that system, which makes calculations simple and transparent. The main motivation of that work 
was to pave the way to a QCD based quantitative study of heavy quarkonium states in the quark gluon plasma 
(several works in this direction have recently appeared in the literature 043); which share with the hydrogen atom 
a number of important features. The main qualitative difference, as far as the bound state dynamics is concerned, 
Oh| between heavy quarkonium states and the hydrogen atom is that vacuum polarization effects in the latter are very 
suppressed. For muonic hydrogen, however, the vacuum polarization effects provide the leading corrections to the 
' Coulomb potential, as it is the case for heavy quarkonium states. This is our main motivation to study muonic 
hydrogen in detail here. 

Muonic hydrogen is under current research at the Paul Scherrer Institute in the so called Muonic Hydrogen Lamb 
| Shift experiment Q- It allows for precision tests of QED which, among other things, probe the electromagnetic 
structure of the proton 13 or the size of the proton Theoretical calculations, on one hand, have achieved an 

impressive precision [lll - llq ] , and a number of experimental results are available [l?], [H| • However thermal effects on 
this atom due to black body radiation, or to electron-positron plasmas, have not been considered to our knowledge, 
either theoretically or experimentally. Current experimental facilities may now produce electron-positron plasmas 
[Tg| . which are also the target of intense theoretical studies [13, (see [Hj for a recent review). Making muonic 
hydrogen atoms slowly travel through an electron-positron plasma would be an ideal experiment to probe how well we 
understand thermal effects in non-relativistic bound states at relatively high temperatures. Recall that an analogous 
experiment at lower temperature with blackbody radiation on Rydberg atoms [23j | first detected thermal level shifts, 
in the early eighties. This would mean a further step in taking advantage of the similarities of the electron-positron 
plasma with the quark-gluon plasma (see [H], [H| for reviews) in order to learn about the non-trivial properties of the 
latter, as it has already been advocated by some authors [26| . 

In the center of mass frame, the proton of a muonic hydrogen is essentially at rest and the muon moves at small 
velocities v < a -C 1. Hence, the relevant scales at zero temperature are those of a non-relativistic system [27|: the 
muon mass m M (hard), the typical momentum p which is of the order m^a/n (soft) and the energy E ~ m M ev 2 /n 2 
(ultrasoft), n being the principal quantum number. Unlike the hydrogen atom case, vacuum polarization effects 
introduce a new scale in the muonic hydrogen atom, the electron mass m e , which is of the order of the soft scale for 
the lower lying states (n = 1,2), but larger for the remaining ones. At finite temperature further scales are introduced, 
not only T, the temperature, but also eT ~ mo, the Debye mass, and others that will be discussed later. In order 
to efficiently deal with the physics at each of these scales we will use the effective theories of non-relativistic QED 
(NRQED) [23|, suitable for energies much smaller than the hard scale, potential NRQED (pNRQED) [28], suitable for 
energies much smaller than the soft scale, and Hard Thermal Loop Effective Theory (HTL) [!§] , suitable for energies 
much smaller than the temperature, in a way analogous to Ref.Qj. Recall that pNRQED facilitates enormously the 
iteration of the Coulomb potential in ultrasoft contribution (at the scale E) , and the HTL action does the same for 
soft thermal photon resummations at the scale eT. When the scales eT and E coincide the combination of pNRQED 
and HTL is crucial in order to obtain consistently both the iteration of the Coulomb potential and the resummation 
of soft thermal photons. 

We will use the real-time formalism [3(| , which is mandatory for the study of the propagation of a (non-thermalised) 
non-relativistic system in a thermal bath. We shall restrict ourselves to temperatures much smaller than the muon 
mass, and hence the thermal bath does not affect the free muon propagator which remains the same as at zero 
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temperature. The same holds true for the free proton propagator, which will be further approximated by that of a 
static source. Thermal propagators will in general be necessary for the photons, electrons and positrons. Recall that in 
the real-time formalism a doubling of degrees of freedom is required to properly account for the thermal propagation. 
External propagators can only correspond to type "1" fields, vertices contain either type "1" fields or type "2" fields. 
Propagators can be "11", "12", "21" or "22". When drawing Feynman diagrams we will understand that all possible 
types of vertices and propagators compatible with a given diagram are added up, and will not display each type 
explicitly (for muons and protons only the "11" propagator must be considered). The techniques and results we shall 
use have been reviewed in ref. [3l| . We reproduced the basic ones in the Appendix [Bl 

We distribute the paper as follows. In Section II we study the ideal case in which the electron mass is set to zero 
(m e = 0). This not only makes calculations simpler, but also makes the system closer to the heavy quarkonium case. 
In section III we focus on the actual case m e ^ 0. These two sections are divided in subsections in which the cases 
T <C p, T ~ p and T 3> p are studied, p being the typical relative momentum in the bound state (or the inverse Bohr 
radius). Section IV is devoted to discussion and conclusions. 



II. m e = CASE 



Let us first consider an ideal case in which the electron mass m e is taken to be zero. This case is in fact closer to the 
one in heavy quarkonium states, particularly in charmonium, than the actual case with m e ^ 0, which we will study 
in the next section. It has already been discussed in the past in order to clarify subtle issues on the renormalization 
group structure of non-relativistic effective theories 32[. 



A. p>T 



For temperatures much smaller than the soft scale (m M a in this case) we can study the thermal effects in the atom 
starting from the pNRQED Lagrangian at zero temperature, up to exponentially suppressed contributions ~ e~ p / T . 
This means that the potentials are the same as the ones at zero temperature. The only difference with respect to 
the hydrogen atom case is that these potentials contain now 0(a) corrections due to vacuum polarization effects 
produced by virtual electron-positron pairs. Like in the hydrogen atom, the ultrasoft photons, electrons and positrons 
are responsible for the finite temperature effects. Our starting point in this section is then eq. (6) from (33j . 



L p nrqed= -J d 3 x^F^(t,x)F^(t,x) + J d 3 rd 3 nS^t,r,R)^id + ^- + j^ + 

+ ^ (- t + 4 *) <3(r) + ics ^k '■ (w x v ) ) S{L r ' R) (1) 

+ J d 3 rd 3 RS^(t,r, R)er • E(t, K)S(t, r, R) + J d 3 xe(t,x)i^D t ,e(t,x) . 

where S(t, r, R) is the muon wave-function field, r being its distance to the proton and R the position of the proton; 
e(t, x) the electron Dirac field, a = e 2 /47r is the electromagnetic coupling constant, and cq, c s and di are matchings 
coefficients which can be found at one loop in (34|. 

Let us separate the cases T < E and T >• E, which are analysed the two following sections: 



1. T<E 



In this case, the leading temperature-dependent contributions are given by the diagram in fig. [TJ in an analogous 
way to the hydrogen atom case. Virtual ultrasoft electron-positron pairs give rise to O(a) corrections and no soft 
thermal photon resummation is necessary at the scale E. Hence, there is no qualitative difference with respect to 
the hydrogen atom, and we will not further discuss it. We refer to [l| for the relevant formulas for the spectrum and 
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FIG. 1: Muonic hydrogen atom self-energy at leading order (due to ultrasoft photons). The double line stands for the propagator 
of the atom and the wavy line for the thermal propagator of the transverse photons. The vertex corresponds to the dipole term 
in the multipole expansion, see |T]). 

decay widths . 

2. T > E 

In this case the scale T can be integrated out before calculating the spectrum and decay widths, we call the 
resulting effective theory pN RQEDyx (for further explanations about the notation see Appendix [XJ. In the photon 
and electron-positron sector this gives rise to the HTL action [29|. In the atom sector, the pNRQED Lagrangian gets 
additional temperature-dependent potentials. At leading order (LO) in a , they arise from the diagram in fig. [T]upon 
expanding E — H in the integrals, and have been calculated in [l| 

^ - i£ - (i + * (&) + i + + °(^f), 

H ~ E ~ a/r. Possible 0(ar 4 ET 4 ) terms arising from higher orders in the multipole expansion cancel out. Note 
that the dominant contribution above is a constant mass shift, and the r-dependent part is (m^a 2 /n 2 T) 2 suppressed. 
Hence, vacuum polarization corrections to the photon propagator may compete with the r-dependent part of the LO 
potential displayed above and must be calculated. From the diagrams in fig. [2] we obtain the next-to leading order 
(NLO) in a, 

5V± NLO) = -§f C (3)TmV 2 + lJ ^r 2 (l + 7 + lo g7 r - log £ + § - 41og2 - 2§|) + (3) 

"l$e+c + 27 + 21og(#))+0 

with m 2 D = (eT) 2 /3 and c a numerical constant. The computation have been done in dimensional regularization 
with e = (4 — d)/2 — > 0. The first line of this result also appears in an analogous calculation that has already been 
carried out in the static limit of the QCD case @. The second line is subleading. We have displayed it to match 
the precision of (jSJ when eT ~ E. In order to eventually check the correct cancellation of the 1/e poles we have 
calculated analytically the leading IR behavior in the appendix [Cj The constant c remains unknown. We see that 

indeed SV^ NL °^ competes in size with SV^ LO \ except for the global energy shift given by the first term in ([2]). In fact, 
it provides the dominant term in the potential for the energy of the photon transitions between two states belonging 
to this case. The IR divergences arising above are canceled by UV divergences divergences arising from contributions 
at smaller scales (E, eT, . . .). These contributions are hard to calculate in the general case because HTL propagators 
must be used for the ultrasoft photons and the Coulomb potential must be kept unexpanded in the atom propagator. 
At these scales, however, the Bose distribution can be expanded, which simplifies somewhat the calculations, and 
produces the so called Bose enhancement, see (IB8|) and (|B10[) below. The dominant contribution arises from the 
diagrams of fig. [3j We have only been able to work out an analytic expressions for the cases eT <C E and eT 3> E, 
which will be discussed below, and for its UV behavior. Technical details for the latter are shown in the Appendix [Cl 
Before discussing the two cases which allow to proceed further analytically, we display the energy shift and decay 
width induced by the temperature-dependent potential ([2]) and ([3]) at leading order in perturbation theory, 

SE n = f£ - 3^l^(0)l 2 (Hl°g(^r) + | + l°g(2T)) (4) 



1 In the case E>T the formulas presented in [lj provide the leading contribution for T 3> ctE only. For T < aE additional contributions 
exist 13511 



FIG. 2: Muonic hydrogen atom self-energy at next-to- leading order (due to ultrasoft photons, electrons and positrons). The 
solid line stands for the thermal Dirac propagator of electrons and positrons and the dashed line for the propagator of the 
longitudinal (Ao) photons. The remaining lines are as in fig. [1] 



-|fC(3)T^ (r 2 >Ti + Q ; 

5Tl = _ aT mUr* )n (l +7 + l 0g7r _ l0g g + |_ 4l0g2 _ 2 m ) + (5) 
+ ^(i + C +27 + 21og(^))+o(^,^), 

in which (r 2 ) n — 2w ^ a 'i [5n 2 + 1 — 31(1 + 1)], n and I being here the principal and angular momentum quantum 

numbers. The labels n, m, . .. are also used through the paper as a short hand notation for the whole ensemble of 
quantum numbers of a given Coulomb state, either bound or in the continuum, and </> n (0) is the wave function at the 
origin. The contributions above are to be added to the ones coming from lower scales, which we discuss below in two 
particular cases that share the feature that the lower energy scales are hierarchically ordered, and hence the method 
of the integration-by-regions can be used [3(| 

• E > eT 

In this case the loop integral is dominated by energy and momentum ~ E for which eT can be treated as a 
perturbation. At leading order in the HTL expansion we obtain 

se; . I yj M , M - Em) (i + ,„ g (^) + 1 - 7 + logW )-^>„ +f , (=£4) 

(6) 

„ p 4a 3 T aTm 2 ^,, , ll9 /l \E n - E m \ 11 , A / ar 2 Tm%\ , , 

«f -^ + ^ JZ ElHr|m)| 2 f - - 2 log 1 " m| +y -log4- 7 + log(7r)j +o( ~ 2 D j . (7) 

We observe that the leading order infrared divergences appearing at the scale T in (|4]) and ([5]) are canceled by 
the ultraviolet divergences in © and (O respectively. Note that the subleading infrared divergence at the scale 
T, in the second line of (O, is very much suppressed in this case (m^a 5 3> am 2 J /m fJ _). A similar calculation for 
the heavy quarkonium case has been presented in [38(. Technical details can be found in the appendix (jC 2[) . 
We only mention here that a collinear region exists that contributes at this order. 

Upon summing up the contributions from both energy regions, namely Q and ^ for the energy and ([5]) and 
for the decay width, finite results are obtained at the desired order, 

SEn = ^ + If \M0)\ 2 (log (flj) - 7 ) + If £ m \(nH m )\ 2 (E n E m )\o gj7 ^ (8) 

cTml(r 2 ) n ^ 3C(3) + ttA + ^ ^ o^Tm*, \ 

JT„ = ^ + (log - 7 + § + log 2 + $f ) _ ^ Em I (n|r | m) |2 log JS^JaL (9) 



FIG. 3: Further contributions to the muonic hydrogen atom self-energy when p 2> T 2> E. The wavy line and the dashed line 
with a blob are the HTL propagators for the transverse and longitudinal photons respectively. 



• eT 3> E, In this case the loop integral is dominated by energy and momentum ~ eT for which E can be treated 
as a perturbation. At leading order in the energy expansion we obtain a contribution which is equivalent to 

-.2 J2 



adding a new term to the potential that goes like m D r 



8E? L ° = ^<r 2 )« + 0(a^mlT), (10) 



STf LO = ^Tn^ {r2)n{ l_ 7 + log?r + bg l_ + 5 } + 0{a 2 r 2 mlT) , , , , 

O 6 TTl O 



We observe that the leading order infrared divergence appearing at the scale T is cancelled by the ultraviolet 
divergence above. An analogous contribution has also been calculated in the static limit of QCD [6]. 

At next-to-leading order in the energy expansion we have restricted ourselves to compute the ultraviolet diver- 
gence analytically, 

^o_^ + c ,_ 2]og ^ +om (12) 

It cancels the infrared divergence in the second line of ([3]), as it should (c* is an unknown constant that can be 
of order 1/a 1 / 2 because of Bose-enhancement). 

Fortunately the contribution of the loop integral for energy and momenta ~ E is subheading. The calculation 
at that scale may even require non-perturbative techniques if E gets close to the scale e 2 T [30l . l39l l40j ] . 

Summing up the contributions from the T energy region and from the energy region, namely (j4]) and (|10p 
for the energy and §5§ and (ITT1) for the decay width, the leading thermal effects for this situation are obtained, 

SE n = ^ + ^M^ + 0(ar^T) (13) 



T 1 . , C'(2)\ 2am 2 n (, T 



(14) 



B. T ~p 



Since m M 3> T still holds, our starting point is the NRQED Lagrangian at T = [27[ (this is correct up to 
exponentially small contributions ~ e~ m >*/ T ). 

£ = ^(^ + ^ + c F e^- + c D e^4 + (15) 

+?cse ^(DxE-ExD) + n+w0n _ i + d* FD 2 F »u + B , D 

8m;. 4 raf, 



FIG. 4: Leading correction to the Coulomb potential due to vacuum polarization when T < p. The thick and the extra-thick 
solid lines stand for the non-relativistic propagator of the muon and the static propagator of the proton respectively. The 
remaining lines are as in fig. [2] There is also a leading order correction coming from an extra diagram, which is obtained from 
the second one by changing the muon line by a proton line. 

where ip and N are the muon and proton Pauli-spinor fields respectively. 

Since T ~ p ^> E, we can integrate out the scales T and p, which leads to what we call pN RQEDt, the suitable 
effective theory for the scale E, which is similar to the one that has already been introduced in section II. A, with 
the only difference that now T is of the same order as the cut off of the effective theory. In the photon and electron- 
positron sectors we have the standard HTL. In the atom sector temperature-dependent potentials are induced. Recall 
that at the scale T there is no enhancement and vacuum polarization effects due to electron-positron pairs are always 
suppressed by a. Hence, the leading potential will still be the Coulomb potential but the first a correction to it will 
already be temperature-dependent. This is given by the diagram in fig. [3J 

The temperature-dependent part of the potential induced by the diagram in fig. [4] is both UV and IR safe in 
momentum space. However, when it is Fourier transformed to coordinate space an IR divergence is encountered. The 
calculation is rather involved so we only display here the final result, which can be given in terms of one-parameter 
integrals of special functions. Details are given in the Appendix [D] (the formulas in that appendix have to be used 
setting m e = for this case). We obtain, 

SVr = - gC(3)Tm^ 2 + (16) 

+^rk Io°° T^TT) M - W + (A 2 + 4) cos(H + 
+pu sin(jdtt) + (6pu + p 3 u 3 )Si(pu)) + 

+ + 7 + log 7T + 2l0gM - 1) + 

_^(i_l 0g(rT) _l 0g7r) + 
+fS££ / °° W rin(pu)(Li 2 (-e«) + «log(l + e«) + £ - £ ) + 0(a 3 T), 
where p = 2rT and Si stands for Sine Integral 

z sint , . 

— dt. (17) 

The LO energy correction is obtained by computing the expectation value of the potential of (fT6]l for the desired 
state and adding the ultrasoft contribution. The calculation in pN RQEDt is identical to one carried out in the 
second case of sect. A (since E ~ m^a 2 , T ~ p ~ m^a, we have eT ^> E). Hence the outcome can be directly read off 
equations (|10[) and (1111) . Notice that the infrared divergence in (|16l) at first order in quantum mechanic perturbation 
theory, induces a contribution that cancels out the ultraviolet divergence in (jlll) . as it should. Also, the contribution 
from integrating out the scale m d can be encoded in a correction to the potential, summing it up to 8V r the following 
finite result is obtained, 

6V = _^_ 3^(3)7^2 + £E&rl + (18) 

+ Jo°° W^TT) (" 4 " 4 A 2 + (A 2 + 4) cos(H + 

+pu sin(yf?u) + (6pu + p 3 u 3 )Si(pu)) + 

§ (- hg(rm D ) + | - 7) + 




FIG. 5: Leading correction to the Coulomb potential due to vacuum polarization when T>p. The lines are as in figs. [4] and 

El 



-^#(i-log(rT)-lo g7 r) + 
Jo" W suu»(Lz 2 (- e «) + ulog(l + e«) + £ - + 0(a 3 T), 
and hence, at first order in perturbation theory, 

8E n = (n\mV\n) (19) 
ST n = -2(n\$8V\ri) 

C. T>p 

Since m p 3> T still holds, we can also start from NRQED at T — 0. Now we may proceed by sequentially integrating 
out first the scale T and next the scale p. After integration of the scale T we get an effective theory which consist of 
HTL contributions in the photon and electron-positron sector, and of NRQED with temperature-dependent matching 
coefficients in the atom sector. This N RQED>t in the atom sector is identical to the one that we have in the 
hydrogen atom case [l| , up to order a corrections induced by the electron-positron vacuum polarization. 

The next step is to integrate out the energy scale p, namely to match N RQED>t to what will be called 
pN RQEDt^T, which is expected to produce temperature-dependent potentials. These potentials must be calculated 
using HTL photon propagators. Let us separate the two following cases: 

• eT~p, 

In this case the computations can be carried out as in section V.B of The relevant diagram is similar to fig. HI 
but instead of a photon propagator with a self-energy insertion we would have to use the tree level HTL photon 
propagator (fig. [5]). The only difference with respect to [l[ is due to the fact that the electron-positron pairs 

that generated the HTL photon propagators are now taken to be massless, so, in fact, the outcome is simpler: 

2 

the non-trivial function g(rn e /3) reduces to 16a j? 2 . We obtain then the following leading order potential 

-m D r f aT 2 \ 
V(r, T) = am D + iaT<i>(m D r) + O , (20) 



where 

dzz 



<t>{x) = 2 



(z 2 + l^ 2 



sin(za:) 



(21) 



This potential coincides with the one first obtained in [2| for QCD (up to trivial changes in color factors made 
explicit in 1]). As in the hydrogen atom case, we can use this result in order to estimate the dissociation 
temperature, which is Td ~ m^a 2 / 3 / In 1 / 3 a, as anticipated in [l| for the QCD case. 

• eT > p, 

In this case T 3> Td always holds. Hence, the imaginary part of the potential is bigger than the real one, so it 
does not make much sense to speak about bound states anymore. 
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III. m e + CASE 

We address now the actual case of a non- vanishing electron mass. Although the real muonic hydrogen is not as 
close to heavy quarkonium systems as the ideal one (m e — 0) , it may still be useful to learn about certain aspects of 
it. In particular about the role of a finite charm mass in the bottomonium system, which is analogous to that of a 
finite electron mass in muonic hydrogen This case may then shed light on the effects of the charm quark mass 
in bottomonium at finite temperature, specially when charm quarks are thermalized. Irrespective of that, muonic 
hydrogen is a real system thai appears in nature, which nowadays is produced in large samples ft and, therefore, 
our results may eventually be checked against experiment. 

For actual muonic hydrogen m P , ~ p for the lower lying states (n = 1,2), whereas for the remaining states (n > 3) 
one may safely consider m e » p . 41| . Let us then analyse these two cases separately. 

A. Lower lying states (n — 1, 2) 

As mention above, these two states fulfill p ~ m e , and hence relativistic electron-positron pairs must be integrated 
at the same time as the momentum transfer p is. Let us see in the following sections how this is carried out depending 
on what the temperature is. 

1. p>T 

Like in the massless case our starting point can be pNRQED. However, now, due to the fact that m e ~ p (rather 
than m e = 0), the electron-positron pairs have already been integrated out when calculating the potentials, and hence 
are not active anymore. The situation is then totally analogous to the hydrogen atom, the only difference being that 
the potentials get 0(a) corrections due to virtual electron-positron pairs, the most important of which is the Uehling 
potential. In other words, the thermal bath contains neither electrons nor positrons, so the thermal effects are only 
due to the photons which do not distinguish between electrons and muons. Hence the results concerning this case can 
be read off section III of ref. [l| by making m — > m M [up to 0(a) corrections]. 

2. T ~p 

Again, like in the massless case, we can start with NRQED. The scales T and m e must integrated out at the same 
time as the energy scale p. In the photon sector, which is not sensitive to the scale p, we get the mass dependent 
HTL action (see section V.A.I of p|). In the electron and positron sectors, which are not sensitive to the scale p 
either, we get a NRQEDt Lagrangian for each of these particles (see section V.A.2 of [l[). In the atom sector, the 
potentials depend now on both temperature and the electron mass, except for the leading Coulomb potential. The 
most important correction is a kind of temperature-dependent Uehling potential, which is obtained from the diagram 
in fig. H 

s = _ 4aV(m.g)m;r _ 2^ rao du ,j _ CQs{(Tu) _ auSl{(Ju)) + (22 ) 

+ <si f 00 duV^+l ( 2 - 3a 2 u 2 + (ct 2 u 2 - 2) cos(™) + 

+ausm(au) + a 3 u 3 Si{au)) - ^Tm 2 D r 2 (l3m e ) 3 J™ ^ u ^^ + 



■ 2rr3 _s 

3tt 



+ 7 + log* + logM 2 - 1) - 



-^TT) (l " ^(rT) - log 2 - (e*- + 1) / °° + ST ^ 

+ / °° $ sm(au) (Li 2 (-e^^+ I ) + (/3m e )V^+l log(l + e^'^+ I ) + 



+ 



+ —^—J0 IF \ „g meV /„2 +1 7 1 e f-c+l + 37r(e^e +1)1,7 l+7+^Og(r/X) +10g7T 



'■+1+1 
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where 8 = 1/T, a — 2m e r and 



(23) 
(24) 



f(m e 8) = / ^ ; — 

r°° x 
q(m e 8) — 8 m„ / da; , . 

Further expressions for these functions can be found in the Appendix B of [l]. The computations that lead to 
are carried out in Appendix [Dj 

Notice that (I22[) has two infrared divergences, which, as in the m e = case, arise when the Fourier transform of the 
momentum space potential is taken, in order to get the coordinate space potential. The IR divergence in the fourth 
line of (|22]l is similar to the one that appears in equation (fT6| for the massless case (with T 2 g{m e 8) instead of m 2 D ). 

2 

The IR divergence in the second last line of (|2"2"|) . however, is proportional to e/j ™* +1 , and hence distinct of the m e = 
case. It emerges from a region in which not only the three-momentum transfer is small but also the component of the 
three-momentum of the electron-positron pair in the loop parallel to the momentum transfer is small. In either case, 
these IR divergences should cancel against UV divergences in the pNRQED calculation. 

The relevant diagram in the pNRQED calculation is again fig. [3J in which the photon line must be understood 
as the mass-dependent HTL propagator (see Appendix [E}. In the dominant contribution to this diagram, E — H in 
the atom propagator can be treated as a perturbation (recall that E — H ~ m^a 2 <C eT). Then using the fact that 
An(fc°,k) is symmetric with respect to fc° — > — k , we obtain 



SE1 



-(r 



0{ar 2 Em 2 D ), 



ST 



eT 



16|_3L 5 (m e /3)(i - 7 + logvr + log ^ + f )+ 



(25) 
(26) 



+ 3,(^+1) (7 - 21o § if + f + lo g( 4 -) - 7 - 21og2) 



0(ar 2 Emj 



We see that indeed the UV divergences above cancel those of ((22]) , as expected. There is a subtle point in this 
calculation, however, which we discuss in the Appendix [El that must be correctly dealt with in order to get the UV 
divergence of the last line (that cancels the IR divergence in the second last line of (|22jl). As in the massless case, the 
contribution from the scale mo can be encoded in a correction to the potential, and this summed to SV ri 



SV = - 



4a /(m e /3)m e r 



2a 2 r« 
irr JO 



(1 — cos(ctu) — auSi(au)) + 



«i roc ^v^+T (2 _ 3a 2 u 2 + {a 2 u 2 _ 2) cos{au) + 



+ausm(cru) + a 3 u 3 Si{au)) - ^Tm 2 D r 2 {Bm e f J Q C 



3 f°° duuVu' 2 + l 



il6a 2 T 3 g(m c P)r 2 



3tt 



iia 2 T 



log(rT) - log 2 - (e 



■(-log(rm D ) 



-7)- 



(in 



due 



+ ^ r °° ^ sin(CTu) (Li 2 (-e^ m « v " 5 + T ) + /3m e V^TTlog(l + e ^ m -^+ I )+ 
_ (^ (u a + i) _ 5(m ^) + i^l) + 

, z4a 2 T f°° <2u [ 5me(g-u)--l _ Sinc(cju) 
+ i JO 7 



+ 1 



28a m~Tr / i / \ 1 4 \ 

fa(e^Wl) I" lo §( r "^ ) + 3 - 7) 



■ ^ ;i) r(-2/3) + 0(a^. 



(27) 



The energy shift and the decay width at first order in perturbation theory can be obtained from ([19 



5. T »p 



This case is very similar to the m e = one. We start with NRQED and integrate out the scale T first. Since 
T 3> m e , the mass-dependent HTL propagators may be expanded in m e /T, and hence become the usual HTL 
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propagators with next-to-leading order contributions due to the non-vanishing electron mass . Hence the finite mass 
effects do not affect the gross features of the system. In particular, in the eT ~ p case the dissociation temperature 
will be similar to the one in the massless case. For eT ^> p, like in the massless case, no bound state is expected to 
survive. 



B. Higher energy states (n > 3) 

As mentioned before, these states fulfill m e » p » Ji, and hence relativistic electron-positron pairs may be 
integrated out before the momentum transfer p is. Let us see in the following sections how this is carried out 
depending on what the temperature is. 

1. m e >T 

In this case we can start with a NRQED Lagrangian for the muon in which the relativistic electron-positron pairs 
have already been integrated out, which gives rise to 1/m^ corrections to the Maxwell Lagrangian. Since there are 
neither electrons nor positrons in the thermal bath, the different situations coincide with those of the hydrogen atom, 
and the relevant expressions can be read off sections III and IV of [l| by making m — > m M (up to 1/ m 2 corrections) . 

2. m e <T 

In this case we can start with a NRQED Lagrangian for the muon keeping the relativistic Dirac Lagrangian for 
the electron in it. When m e ~ T both scales must be integrated out at the same time. In the photon sector, 
a mass-dependent HTL Lagrangian is induced (see section V.A.I of 1]). In the electron and positron sectors a 
temperature-dependent NRQEDx Lagrangian for each particle is induced (see section V.A.I of [1]). In the muon 
sector, a temperature-dependent NRQED>y Lagrangian is also induced. At lower orders, it can be obtained from the 
diagrams of section IV of [l[, together with diagrams containing an electron-positron loop. The most important effect 
is the appearance of two mass shifts, one ~ aT 2 /m^ from the diagram (38) of [l| and the other one ~ a 2 m e from 
the second diagram in fig. [4] In the proton sector an analogous mass shift ~ a 2 m e occurs, that is due to a diagram 
obtained from the previous one by changing the muon line by a proton line. 

The next step is to integrate out the scale p, namely matching to pNRQED using the HTL Lagrangian. This has 
already been done in section V.B of [l| - At leading order, this produces a temperature-dependent potential and further 
mass shifts, which can be read of a corrected version of formulas (58) and (60) in that reference (making m — > m e ). 
The origin of the corrections is discussed in Appendix [El and boils down to a simple replacement, see (|E17I) . Both the 
potential and the mass shift contain an imaginary part, which becomes more important than the real part starting at 
some temperature T^, which we call dissociation temperature. In the following section the dissociation temperatures 
will be estimate for several states. 

For T m e none of the higher energy states exists anymore, since the dissociation temperatures fulfill m e > T, 
see next section. 

C. Dissociation temperatures, level shifts and decay widths 

The dissociation temperatures will be estimated in a similar way as they were in the hydrogen atom case [l[ . This 
is by identifying the momentum scale for which the real and imaginary part of the momentum space potential are 
equal, p ~ (16a) l / 3 (g(m e 0) + (m e {3) 2 riF(m e l3)/2) 1 / 3 T =: md, and equating it to the typical momentum transfer in 
the muonic hydrogen atom 2 , p ~ m^a/n 2 . The results are displayed in table I. Table II shows the same results for 
the hydrogen atom. In order to carry out the estimates we have use the formulas of Appendix [El which are meant 



2 Note that in ref. [J the typical momentum of the electron, p ~ m e ct/n, was used rather than the typical momentum transfer, 
p ~ m^a/n 2 , as we have adopted here. For the lower lying states the order of magnitude estimates does not differ much, but for higher 
energy ones may differ considerably. We have also identified an error, that affects the dissociation temperatures, in the computation of 
the potential (57) of [lj, which is explained and corrected in Appendix Q5] We reproduce here the table I of ref. [l[ in table II with the 
current choice of p and the correct version of the potential for the sake of comparison. 
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E (MeV) 
-0 . 0028 



-0 . 0032 
-0 . 0034 
-0 . 0036 
-0 . 0038 



FIG. 6: E vs T for the IS state. This result was computed using the assumption that - 3> ran, so we expect important 
deviations from the real energy starting at T ~ 4MeV 

for higher excited states (T ~ m e ^> p), for the case of the lower lying states, instead of the m e = formulas. This is 
indeed legitimated: We are just not taking advantage of the fact that T<* 3> m e for these states, which produces the 
simplifications discussed in section UlI A 31 
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T d (MeV) 


m D (MeV) 


m d (MeV) 


1 


1.7 


0.16 


0.77 


2 


0.41 


0.036 


0.19 


3 


0.19 


0.012 


0.086 


1 


0.13 


0.0056 


0.048 


5 


0.10 


0.0030 


0.031 



TABLE I: Dissociation temperature for the lower lying states of muonic hydrogen 



11 


T d (KeV) 


m D (KeV) 


m d (KeV) 


1 


49 


0.15 


3.7 


2 


36 


0.020 


0.93 


3 


31 


0.0061 


0.41 


4 


28 


0.0025 


0.23 


5 


26 


0.0013 


0.15 



TABLE II: Dissociation temperature for the lower lying states of hydrogen 

From the results in table I, we see that only the lowest lying states n = 1, 2 survive at temperatures of the electron- 
positron plasma m e < T Hence, only transitions between these levels might be observed in an eventual experiment. 
We shall focus on the experimentally prominent K a transition [T7l |. We display our results for the energies of the 15 
and 2P states, and for the energy of the K a transition as a function of temperature in the range T £ (2, 0.05)MeV. in 
fig. [51 fig- [7]and fig. [^respectively. The calculations have been carried out numerically, using first order perturbation 
theory for the potential (|27p . 

Another observable that we can predict with our results is the decay width. In fig. [9] we show the decay width for 
the IS state as a function of the temperature. Comparing fig. [5] and fig. H]it can be seen that Td ~ 1.7 'MeV is the 
temperature that makes the decay width of the same magnitude as the binding energy. 
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E (MeV) 
-0 . 0007 

-0 . 0008 

-0 . 0009 



-0 . 0011 



FIG. 7: E vs T for the 2P state. For the 2P state the typical radius is four times bigger than for the IS state, so we expect 
important deviations from the real energy starting at T ~ lMeV 



E (MeV) 
.0025 

. 0024 

. 0023 

. 0022 

. 0021 



0.0019 



FIG. 8: K a transition vs T. This result is reliable until T ~ lMeV. Note that this temperature is twice what we predicted for 
the dissociation temperature of the 2P state 



Decay Width (MeV) 






FIG. 9: Decay width for the IS state. As for the binding energy this has been computed using the assumption that - S> hid 
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IV. DISCUSSION AND CONCLUSIONS 

In this paper we have discussed the properties of muonic hydrogen in a thermal bath, which may consist not 
only of blackbody radiation but also of an electron-positron plasma. We have further developed the effective theory 
techniques for bound state systems at finite temperature initiated in [![, in particular the application of dimensional 
regularization to the factorization of the various scales in the system. They facilitate enormously the organisation 
of the calculation. For instance, they make apparent when Coulomb or HTL resummations are necessary and when 
they are not. In addition, both partial and final results are naturally obtained as a series of small scales over large 
ones, thus providing a good control on the systematics. 

We have discussed two cases. We have first addressed the academic case of muonic hydrogen with a vanishing 
electron mass, which turns out to be closer to heavy quarkonium states than the actual case with a non-vanishing 
electron mass, that we have addressed next. All the thermal modifications we have found turn out to be spin 
independent. 

In the zero electron mass case, we have studied how the effects of vacuum polarization modify the picture that we 
encounter in normal hydrogen [1[. The modifications turn out to be important when the temperature is larger than 
the binding energy. For instance, they would give the leading order contribution to a hypothetical K a transition for 
high enough temperatures ©■ For temperatures below dissociation, we have presented the leading order, and selected 
next-to-leading order, thermal corrections to the binding energy and decay width. 

In the actual electron mass case, muonic hydrogen behaves very much the same as hydrogen for temperatures below 
the electron mass. For temperatures larger or of the order of the electron mass the vacuum polarization effects are 
sizable, and, at some point, make the bound states dissociate. We display in table I the dissociation temperature 
for the lower laying states. We have also calculated the thermal modifications to a number of observables before 
dissociation occurred. For instance, we plot the dependence of the K a transition on temperature in fig. [31 which 
could be tested experimentally in the future [8] . 

We close with a concrete application to the heavy quarkonium case. As we have mentioned before, the way a finite 
electron mass affects muonic hydrogen is similar to the way a finite charm quark mass affects bottomonium [4l| . 
Since this should also be the case at finite temperature, we can easily translate to the QCD case the results for the 
dissociation temperature of muonic hydrogen, which we show in table III. 



m c (MeV) 


T d (MeV) 


oo 


480 


5000 


480 


2500 


460 


1200 


440 





420 



TABLE III: Dissociation temperature for Upsilon (IS) for different values of the charm mass. The nr = 3 light quark masses 
are set to zero. We use as an input the values of the Bohr radius and Aqcd found in table 2.1 of [L|. The values of these 
parameters for uf — 3 are used for all values of m c except for m c = 0, where we use the ones for uf = 4 
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Appendix A: Notation for the different effective field theories 

At zero temperature there are three different energy scales for non-relativistic bound states. These are the hard 
scale (for muonic hydrogen m M ), the soft scale m^a and the ultrasoft scale m^a 2 . Moreover, a finite temperature 
system also has a different energy scale as T, eT, e 2 T... This makes it hard to find a comprehensible notation for 
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all the effective field theories that may arise from integrating out the different degrees of freedom. In this paper we 
have used the following notation. Basically we name the effective field theories as one would do for zero temperature, 
and we encode the temperature information in a subindex. The subindex T means that the temperature has been 
integrated out, the subindex m B means that also the scale eT has been integrated out. Since the matching coefficients 
of the effective field theory that we obtain after integrating out m M (p) and T is not the same if (p) ~ T or if 
to a» (p) » T we include a symbol <, > or blank depending of the relation between these scales. For example, if we are 
in rrif, (p) ~ T we will arrive to NRQED T (pNRQED T ), but if (p) » T we reach NRQED >T (pNRQED >T ) 
because T is smaller than the energy cutoff of NRQED (pNRQED). 

Appendix B: Basic formulas 

In this appendix we display a number of formulas of the real-time formalism that are relevant for the paper. Our 
notation closely follows ref . [3l[ . Recall that in this formalism there is a doubling of degrees of freedom [3(J ■ Fields are 
labeled as "1" or "2". Fields "1" ("2") only interact with fields "1" ("2") according to (minus) the original interaction 
Lagrangian. Fields "1" may be converted to fields "2" , and vice versa, through propagation so that propagators 
become 2x2 matrices. For instance, for a free scalar field we have 

A(K) = ( **-™*+* ° )-2mS(K 2 -m 2 ) ( " s(fco) ^o) + n B (k )\ (m) 
V K*-m*-ie J \8(k ) + n B {k ) n B (k ) j 

where iA(K) is the Feynman propagator, and n B (ko) the Bose distribution function. 

For the tree-level propagator of the transverse electromagnetic field in the Coulomb gauge, iAij(K), we have 

A ij (K)=(s ij -'0^A(K) (B2) 

The tree-level Ao-propagator matrix in the Coulomb gauge is diagonal, traceless and the "11" component coincides 
with the propagator at zero temperature. 

For the tree-level propagator of a Dirac fermion field, iS(K), we have S(K) = (If + m) A(K), where A(K) follows 
from A(K) by replacing n B by —rip, Hf being the Fermi-Dirac distribution. 

Since we are always in the case >> T all vertices involving muons or protons will be of type "1". However, 
vertices involving photons and electrons can be both type "1" and "2". At the order we are calculating, it turns out 
that we only need propagators of the type "11" for the photons, either at tree level or including one-loop self-energies. 

For computations that require loop corrected propagators (for example fig. [3] and fig. [2J it is convenient to use the 
so called Keldysh representation [31(. In this representation the retarded, advanced and symmetric propagators are 
defined as, 





= A n 


- A12, 


(B3) 


A A 


= A U 


- A 21 , 


(B4) 


A s 


-A n 


+ A 22 , 


(B5) 



Notice from (|B1|) that at tree level only As depends on the temperature. 

In order to calculate loop corrections to the An propagator in an efficient way, we use the following method [31(. 

• We compute the A^ propagator, using the fact that for this propagator the Dyson equation is of the zero 
temperature type (this is not so for As). 

A R = A° R + A B U R A R , (B6) 

where the self-energy II# is 

n R = n u + n 12 , (B7) 

and A^j is the tree- level retarded propagator obtained according to formulas (IBip ~ (|B3[) . 

• The advanced propagator A^4 is the complex conjugate of the retarded one and the symmetric propagator As 
reads, in the bosonic case, 

A S (K) - [1 + 2n B (\k \)]sgn(k )[A R (K) - A A (K)}. (B8) 



15 



• Finally, one obtains 

A u = i(A R + A A + A s ). (B9) 

Note that for \ko\ <C T the size of the symmetric propagator for bosons is larger than what one would expect from 
naive power counting 

A S (K) = %-[A R (K) - A A (K)} (BIO) 
fc 

plus terms suppressed by 1/T. This effect is called Bose-enhancement, and it complicates the power counting of the 
EFTs at scales lower than the temperature. 



Appendix C: Calculations in Section II. A 

All these computations have in common that the starting point is pNRQED. In this effective theory, and in all 
the theories that are derived from it by integrating out further degrees of freedom, the leading correction to the 
Hamiltonian is @, [H| 

5H = ieV / £w E-H-k + ie ^ + ^ Ao °^ (C1) 

with A^ and A 00 being the "11" transverse and A photon propagators respectively. SH is a potential SV (energy- 
independent) only if E is much smaller than the scales inside the A propagators. Depending on the concrete calculation 
that one is doing the explicit expression for A may be different. 



1. Integrating out the T scale 

After integrating out the scale T we will reach pN RQED>t- At leading order the difference between pNRQED 
and pNRQEDyr will be a correction on the potential of the type of (|Clj) where the internal momentum K is of order 
T. 

The LO correction @ is obtained using the tree level photon propagator in (|C1[) . If one uses the Coulomb gauge 
thermal effects only appear in transverse photons. This computation is done in detail for different relations between 
T and E in [l|. 

= *> J ( 4 « - w) **<**>«<w «»> 

The NLO potential comes from including one-loop corrections to the photon propagator (this corresponds to the 
diagram of fig. [3]). This can be done analytically for T » E. The required expression for the photon propagator are 
found in section III.B of 6]. Since T >> E the muonic hydrogen propagator can be expanded 

1 , 1 ( E - H ) i... (C3) 

E — H — ko + ie — ko + it (— fco + it) 2 

Let us compute the first term in this expansion 

5V NL °» = ie V / |^-J_ (fcgA^ + fc^-Ajo loop y (04) 



Looking at the expressions of [6( and using the relation (|B9j) one sees that both Ay and Aoo are even functions in fco, 
so one can use for the first term of the expansion 

n8(-k ). (C5) 



-fco + it 

So we only need the photon propagator in the limit fco — > to get the leading order 
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In [y] the propagators is given as an integral of a parameter qo. The best strategy to perform the calculation is to 
integrate first over internal momentum and leave the integration of this parameter to the end. This computation was 
done in section IV.B.l of [6j, here we take the Abelian limit making Cf = 1 and Ca = 0. 

rNLOa _ 3 /-/ Q \ a „2t, 2 , -"2^2 ( 1 , , l l„„T 2 , 2 ,fl„„o r,C'( 2 ) 



^™ iUa = ^( 3 )- r Tm D + i-Q^Trnj} I - + 7 B + logvr - log + - - 4 log 2 - 2^yj (C6) 

For the next-to- leading order in ([3]), namely the contributions coming from second term in (|C3[) . we will restrict 
ourself to the extraction of the leading IR logarithmic behavior. We use the fact that in the infrared the photon 
self-energy approaches the HTL limit, so instead, we substitute the complete self-energy n HT Le~ /3 k . The factor 
e -0 k j g introduced to regulate UV divergences and does not affect the leading IR behavior we are interested in. 
The HTL self-energy have the property that can be written as Hhtl = mr)f(ko/k) with f(x) a non trivial function. 
Performing the change of variables fco = kx the infrared behavior can be easily extracted. First we compute the 
contribution from the retarded part of the longitudinal photon propagator. 

NLOb_-2 r i (J? _ m f d D k fc^ nj^fc) ,_ e 2 f d D k n£(fc ,fc) 



5V L -ier\E H) J — _^-_L_^ _ ___ J ( C 7) 

We have used that r l {E — H)ri = P ms terms that vanish on the physical state. As we are only interested in the 

leading logarithm behavior, the change n^(fc , k) — > n^ Ti R (k , k)e~^ k2 — m 2 D f R (k /k)e~^ 2k2 can be made, where 



= a log ( (cs) 



gyNLOb 



D 



2m,, 



d D k fk(ko/k)e-? 2k2 
{2tt) d k 2 (k -ie) 2 



(C9) 



The equality above is only true as far as the leading IR behavior is concerned. Now with the change ko = kx 



gyNLOb 



ie 2 m 2 D &d-i 
' 2to m {2ir) D J 



dkk D - 5 e-? 2k2 



dxfk(x) 
(x — it) 2 



ie 2 m 2 D £Id-i 



4m M (2tt) d \fi 



T 



D-4 



r(- 



,D 



(x — ie) 2 



. (C10) 

£Id-i is the D — 1 dimensional solid angle. Note that the contribution from the advanced part of the longitudinal 
photon propagator is obtained by replacing fj^(x) — > in (|C10|) . Since it has all the singularities in the upper 

half complex plain the corresponding integral gives zero. The integral over x in (jClOj) can be done using standard 
techniques from integration in the complex plane, and the result is 



dxf^jx) 



(Cll) 



The contribution from the symmetric part of the longitudinal photon propagator does not produce IR divergences in 
dimensional regularization. So we do not need to compute it here as we are only interested in the logarithms. 

Next we proceed analogously for the transverse photon propagator. As in the longitudinal part, only the retarded 
plus advanced contribution is infrared sensible. We approximate n|!(fco, k) — m 2 D f^{ko/k)e~^ k where now, 



, o „ X ( X + 1 



ic 



(C12) 



then 



5V" LOb 



1 {E-H)i 



d D k 



(27T 



j.mle-^ 2 



kikj 



fl(k /k) 



(k 2 — k 2 + isgn(ko)e) 2 {k 2 



1 — isgn(ko)e) 2 

(C13) 

The first (second) term in the square brackets have all the singularities in lower (upper) complex fco half plane, and 
hence the whole expression vanishes. The contribution from the symmetric piece of the propagator also vanishes for 
the same reason. Hence, 



SV 



NLO 



= 



(C14) 
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2. The E scale for T » E » m D 



In order to obtain ([6]) and ([7]) the starting point is pN RQEDt- In this effective theory the atom self-energy gives 
the correction to the Hamiltonian (jCll) . Recall that the photon propaga tors in (|C1|) must be taken in the HTL 
approximation. We use the method of the integration by regions [36l . l37j in order to evaluate it. In this integral 
there are three relevant regions: k, ko ~ E with A = ko — k ~ E (we call it off-shell region), k, ko ~ E with A ~ mD 
(collinear region), and k, k ~ (m£)-region). By power counting one can see that the m^j-region will be just a 
higher order correction, so we focus in the regions k, ko ~ E. In the off-shell region the atom propagator can not be 
expanded, but the HTL photon propagator can. For example, the longitudinal photon retarded propagator can be 
written as 



and, after the compulsory expansion 



Aoo > R - ki + mifkiko/ky (C15) 



_ 1 m%fk(ko/k) 
Aoo,i?. - p + ... (016) 



The fact that the non-trivial functions appears only in the numerator after the expansion is crucial in order to be able 
to do the integration analytically. The collinear region does not contribute in the part that is related with longitudinal 
photons. 

The transverse photon retarded propagator is 

A V> R = U2 U 2 ^ ™2 fTl/. TT^ , 7" ( C17 ) 



\°V fc 2 j 

kl - k 2 - m%j%(ko/k) + le 



In the off-shell region a expansion like the one for the longitudinal photon propagator has to be made. In the collinear 
region (that has k,k ~ E but fcg — k 2 ~ m D ), the atom propagator can not be expanded also, but the HTL photon 
propagator has to be expanded around the region \ko/k\ ~ 1. As in the previous case, this expansion makes it possible 
to proceed analytically. 

Details of this computation can be found in [38j , 



3. Integrating out mo for mo » E 



After integrating out mo we will arrive to what we will call pN RQED >mD . In the photon sector we will have a 
non trivial action but we will not need it at the level of precision we are working. In the atom sector there appears 
a correction of the potential in the matching between pN RQEDt (or pNRQEDyT) and pNRQED >mn of the form 
of (jCip where the internal momentum is of order k ~ mp. 

Because >> E we can also put E — H = at leading order in the atom propagator as in section [C"T1 Hence, 
we will need the HTL photon propagator in the ko — > limit(i.e. T » k » k ). The HTL photon propagators are 
very well known and can be found in 0, [3(J HJ . In our case the fact that ko — > makes 

klAij(K) = 0, (C18) 



/ 1 nirTm 2 \ 

^w^H^' »(^4f )' (C19) 

Using this in (|Cip the results (fT0|) and (|TT|) are obtained. 

For the next-to-leading order (|12l) the calculation is very similar to the one we have carried out for the T scale. Like 
in that case, we will focus on the logarithmic behavior, now in the UV. As an example, we study again the retarded 
part of the longitudinal photon propagator (note that consistency with eq. (IC14|) ensures that there will not be a 
logarithmic contribution from transverse photons) 



gyNLOm D _ - n 2 



L leT J (2^ (ko-ie) 2 k 2 + m 2 D ft(k /k)' 2mJ (2t:)D (k ~te) 2 k 2 + m 2 D f^k /k) 

(C20) 



d D k hk^E-H) 1 , ie 2 f d D k k 2 1 
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using that scale-less integrals in dimensional regularization are zero 



c v NLOm D = ie 2 m 2 D r d D k fjj(k /k) 1 

2m M J {2ir) D (k -ieyk 2 +myii(k Q /k) 1 ' 



Now is when the change fcg = kx becomes useful 



c v NLOrn D _ ie m D ii D _i roo dxf^(x) roo dkk D ~ 3 _ lrOr>\ 
0V L - 2m, (2tt)B J-oo (x-ie) 2 JO k 2 +rn 2 D f^(x) ~ l^ZZJ 

_ ie 2 m 2 D np-i ^g^^Dn^ ^m D ^ £> ~ 4 J°° fcCfjMl^lf 



4m M (27r) D " L " ! ' L V 2 i \ J-oo (x-ie) 2 

If one is only interested in the logarithmic behavior, one can put D — 3 = 1 in the x integration and use (jCll[) . 
Formula (|12[) is readily obtained from the expression above. 

4. UV behavior for E ~ mzj 

Although in this situation we have not been able to obtain an analytic result for (|Clj) . its UV behavior can be easily 
isolated as follows. 

1 1 , (E-H) ^ , 1 (E-H) 



E — H — ko + ie \E — H — ko + ie —ko+ie (—ko + ie) 2 / —ko + ie (—ko + ie) 2 

The piece in brackets in the right hand side of the equation above leads to a ultraviolet finite expression when 
substituted in (jCip . So the ultraviolet divergences arise from the remaining terms in the right hand side of this 
equation. In fact the computation for these UV divergences is exactly the same as in the case mo >> E, which we 
have carried in section [C"3l 

Appendix D: Calculations in Section III. A 

1. Correction to the Coulomb potential in pN RQEDt 

In this part we deal with the matching procedure that has to be done for T ~ p. In perturbation theory the 
pN RQED potential is related to the Fourier transform of the longitudinal photon propagator in the limit where 
Po — > (for P the external momentum of the propagator). For this temperature range the propagator, which can be 
obtained by the procedure outlined in Appendix [Bl is needed for T ~ m e ~ p >> po- The retarded self-energy reads 

u R (p) = -k r ^ + 4- r dfcfc(2fc2+2T "'- p2/2) log ( ^m) - (di) 

2ie 2 p roo dkk ie 2 p m 2 1 

T P J P/2 pj k 2 +m 2 rrp fl J„2 /4+m 2 • 

e V + 1 e V +1 

For ni e — > 0, this self-energy coincides with the Abelian limit of the one found in Q. 

In order to obtain the potential the first step is to use formulas (|B9[) and (|B8|) to get the corrections to the "11" 
propagator. Then the propagator is related with to potential plus the self-energy by the following formula 

V(r) = -e 2 J (|^(e ipr - l)An(po = 0,p) (D2) 
at leading order this gives the Coulomb potential V(r) = — — . For simplicity we define 

/ (Sj^ etPrAll{P0 ^^ P) {D3) 



V r (r) = -e 
and 



r,„ , 2 I d J lp i A ll (p D = 0,p) (D4) 
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such that 

V(r) = V r (r)-V m 
For the next to leading order we need the corrections to the propagator 

6 A n = Ai + A 2 + A 3 + A 4 + A 5 

where 

2e 2 f°° dkk 2 



Ax 



"V Jo y/tf +ml{eN k2 +™l + 1) 



(D5) 
(D6) 

(D7) 



A 2 = 



2e 2 f°° dkky/k 2 

k 2 P 5 Jo ~Jhf&+r, 



■ loe 



+ 1 



-p + 2k\ 
\p + 2k\ 



2e 2 

ir 2 p 5 J 



dkk 2 ^Jk 2 + m 



2 pi 



dX 



1 



1 



p — 2k\ + ie p — 2k\ — ie 

(D8) 



27^ 



A — e 2 P« 
^ 3 ~~ 27r 2 p a JO 



7^ Jo 



dkk 
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dfcfc 2 r 1 fi\( i i ^ 

,/V^+^I+l) J-l UA ^p-2fcA+ l£ ^ p-2k\-iej 



(D9) 



A 4 



4iTe 
irp 5 



2 poo 



dkk 



p/2 e 



(D10) 



2iTe m 



2^,2 



7rp J 



s /3 v 'p 2 /4+m 2 + j 



(Dll) 



The contribution of (|D7[) to (|D2|) leads to the first term in the first line of (J22J) . In order to calculate the contribution 
(|D8|) to (|D2[) . it is convenient to leave the integration over the internal momentum k to the end. Consider then, 



2 // A / (2tt)^ P 5 



1 



1 



p — 2kX + ie p — 2kX — ie 



and use 



p-2kX + ie 2kX-ie {2kX - ie) 2 {2kX-ief (2kX-ie) 4 {2kX - ie) 4 (p - 2kX + ie) 
in dmU), 

£ f 1 ^\ ( 1 4- i ^ f ^"'p e ipr _)_ £ f 1 ,/\ f i + i A r rf "^ e'P r 

2 J-l " A \j2fcA-ie) 2 T (2fcA+ie) 2 J J (27r) B -! p 4 ~r 2 J-l " A \ (2/cA—ie) 4 ~r (2/cA+ie) 4 ^ J (27T) 13 - 1 p 2 

-sir 1 n>A r d °' 1 p e' pr / i | i A 

2 J-l "J (2tt) d - 1 p \^(2feA-ie) 4 (p-2feA+ie) ' (2fcA+ie) 4 (p-2fcA-ie) J 



(D12) 



(D13) 



(D14) 



Some terms vanish because of the symmetry A — > —A. The integral over p in the first and the second term is 
straight-forward. By using the symmetries in the A and p variables, the third term can be simplified as follows, 



32w' 2 r f-1 ^ A /-oo dp(e lpr e lpr ) ( (2feA-ie) 4 (p-2fcA+ic) + (2fcA-He) 4 (p+2/cA+ie) + — 
■•■ + (2fcA+ie) 4 (p-2fcA-ie) + (2/c A-ie) 4 (p+2fcA-ie) ) - 

At this point the integral over p can be done using standard techniques of complex analysis, 



(D15) 
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(D16) 
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Hence, our final result for (|Df 21) reads, 
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8nr 



dX 
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(2kX + ie) 4 {2kX-ief 
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(2k\ + ie) 4 (2kX-ie) 4 
(D17) 

After performing the integration in A we obtain the second and third line in (|22j) . A3 can be computed in a very 
similar way, and leads to the second term of the first line in (f2"2"j) . Let us next consider A4, 
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(D18) 



We have separated above the pieces that lead to infrared divergences from the ones that do not. The first term inside 
the square brackets gives the fourth line in (f2"2"j) . and the second term together with V m gives the fifth line of (|22|). 
The rest of terms in the square brackets give the sixth and seventh line of (|22[) . A5 can be computed in a very similar 
way, and gives the remaining terms of 



Appendix E: Computation of the HTL retarded self-energy for the longitudinal photon in the m e 7^ case 



There are some subtle points in the computation that leads to (126[) . which do not arise in the massless case and are 
worth elaborating upon. Let us start from the formula (37) of [3l| for the massive case, 



-2ie A 



d D k 



(27T) 
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(q Q k Q + qk + m 2 e )[A s (Q)A R {K) + A A (Q)A S (K)] 



(El) 



Q = P — K . It is customary to make the change K — > —Q in the second term to get a simplified expression that 
reduces to twice the first term. However, the terms proportional to rrig, which do not exist in the massless case, have 
a stronger IR sensibility than the remaining ones. This leads to ill-defined expressions in the HTL approximation 
for T ~ m e . These expressions must be properly defined in order to get consistent results before and after the shift 
K — > — Q has been carried out. 

To see this in detail we work out just a small part of the computation which illustrates the point, namely the part 
of Sn^(P) that comes only from the ml term in the numerator. We call this term n m 2(P), 



n m2 (P) = -2ie 2 m 2 



d D k 



(2tt) d 

One can use the shift K — > — Q to get a simplified expression, 

oof d k 



[A S (Q)A R (K) + A A (Q)A S (K)} 
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By carrying out the HTL expansion and taking into account that we must expand also po ( p S> po in the computation 
of the potential), we obtain, 
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Notice that the last term is ill-defined in the IR (this is apparent if spherical coordinates are used). Let us focus on 
the imaginary part, 
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For illustration purposes, we will make this computation in two ways 
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• Using d k = c?0d_2 d(cos8)k dk, where 9 is the angle between the internal momentum and the external 
one. 



(vtt no- 2 2 2 ".D — 2 

iyll m 2 = — ilie m e pop 



and performing the angular integration one arrives at 

P (2tt)^-i Jo 



dkk D - 2 n F (y/k 2 + m 2 e ) 
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(E6) 
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This expression has an end-point singularity, but one can skip it in dimensional regularization by choosing 
D — 4 > 0. So the result is zero in this way. 

• Using d^k = d°- 2 k ± dk z 
We choose z to be the direction parallel to the external momenta. 

O f°° i f c 

Sh m , = -16ie 2 m 2 ePo p 2 y dfcj.fc£- Vty*! + fc? + m 2) y 
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1 
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Notice that the integrand vanishes for all k z except when k z ~ 0, so one may substitute 

dnp k 2 



n F {\jk\ + k 2 z + m 2 ) -> n F (J~k? x + ml) + 
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This simplifies the computation a lot because complex plane integration techniques can be applied, 

~„ 9 9 /"°° rffc i fc^ 1-3 drip /"°° ,„„ „, 

Sn mr Wm Po ^ / ^ f / dk z 8(2pk z ) E10 

(27T)^ 1 7 ^fc2+™f d£ 7_ QO 

This expression does not have any end-point singularity, and in fact it docs not need a regularization anymore, 
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This is indeed the expected result, that eventually leads to terms contributing to the second line of (121)1) . We arrive 
then at the paradoxical situation in which the final result depends on the precise way DR is implemented, either in 
spherical coordinates or in Cartesian ones. The apparent contradiction is resolved by noticing that DR in spherical 
coordinates does not allow for the shift K —Q. In order to show this is the actual reason for it, let us start now 
from cq. (|E2l) 

3II roi (P) = -4ir 2 e 2 m 2 e J ^^S(K 2 ^m 2 e )S((K~P) 2 ~m 2 J[ S gn(k -p Q )(l-2n F (\k \))- S gn(k Q )(l-2n F (\k Q -p Q \))} . 

(E12) 

The integral over fco is straight-forward, 

QU m j(P) = -2e 2 m 2 e n f ^J +m2 [(n F (^k 2 + m\ - p ) - n F (y/k 2 + m 2 ))5(p 2 - p 2 - 2 Po ^k 2 +m 2 e + 2kp) - ... 

... - (n F (y/k 2 +m 2 e + p Q ) - n F (^/k 2 + m 2 ))8(pl - p 2 + 2p Qy /k 2 +m 2 + 2kp) . (E13) 

So far we have used the complete expression for the self-energy. We apply next HTL expansion and also p 3> po, like 
above, 



Sffl m2 (P) = -2e 2 m 2 e n 
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From this expression, no matter if one uses spherical or Cartesian coordinates, one arrives at, 
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which agrees with (|E1 1|) . If the expansion for small po is not carried out the same result is found with the three 
methods, because then no regularization is needed. 

The conclusion is that when the formula for the retarded self-energy in the massive case i.e. 



K L R (P) = 4e 2 / - J „„„,;, ; ' ^ rs , (E16) 

pi 



(2tt) 3 ^/WT^ e /3v /p +^f + 



Po - + it 



is expanded for p 3> po, DR must be used in Cartesian coordinates in order to properly regulate IR divergences (i.e to 
be consistent with the shift of momenta carried out at some point in order to get the expression above). This point 
was overlooked when calculating the potential in formula (57) of [l| , and terms analogous to (|E14|) were missed. The 
correct formula is obtained by making the following substitution in (57) of 

m 2 B 2 

g{ m p)^g{mP)+ 2{e0m P +i) (E17) 
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